Fine structure of chiral symmetry breaking in the QED3 theory of underdoped high- Tc 

superconductors 



Babak H. Seradjeh and Igor F. Herbut 
Department of Physics, Simon Eraser University, Burnaby, British Columbia, 

Canada V5A 1S6 



(N 
O 
O 
(N 



00 



O 
o 

I 



I 

o 



(N 

O 
(N 
O 



I 

O 

o 



We study the effects of tlie leading irrelevant perturbations 
on chiral symmetry breaking in the effective QED3 theory 
of d-wave superconductor-insulator transition in underdoped 
cuprates. For weak symmetry breaking, the effect of a pertur- 
bation on energies of various insulating states can be classified 
according to its engineering dimension in the maximally sym- 
metric theory. Considering the velocity anisotropy, repulsive 
interactions, and higher order derivatives we show that the 
insulating state with the lowest energy is the spin density 
wave. 



I. INTRODUCTION 

The d-wave superconducting state, besides the famil- 
iar rotational invariance, at low energies possesses an ad- 
ditional "chiral" symmetry for its quasiparticle excita- 
tions 0]. It was recently shown |0| that the quantum 
phase transition from the antiferromagnetic insulator to 
the d-wave superconductor may be understood in terms 
of this hidden symmetry: while present in the supercon- 
ducting state, chiral symmetry is manifestly broken when 
the spin density wave (SDW) order develops. Dynami- 
cal agent that brings about this change of symmetry was 
argued to be provided by the vortex fluctuations of the 
condensate, that can be represented by an effective gauge 
field minimally coupled to gapless quasiparticles of the d- 
wave state @, §. That way the theory of phase fluctuat- 
ing d-wave superconductor can be mapped onto the three 
dimensional quantum electrodynamics {QED3), which 
was proposed as an effective low-energy description of 
the underdoped high-temperature superconductors 0|. 

Chiral symmetry, however, is not an exact symmetry of 
the d-wave superconducting state, but arises only in the 
low-energy limit of the standard quasiparticle action. At 
low energies one is allowed to linearize the spectrum near 
the nodes of the d-wave order parameter and neglect the 
higher order derivatives and the short-range interactions 
between quasiparticles, which are both linearly irrelevant 
by power counting. In this approximation the action be- 
comes symmetric under a global U{2) x U{2) transforma- 
tion, where each U(2) factor acts onto a four-component 
Dirac field that describes the gapless excitations near the 
pair of the diagonally opposed nodes [Q (see Fig. |l|). If 
one would also neglect the intrinsic velocity anisotropy 
of the d-wave superconductor vp 3> wa, with vp being 
the Fermi velocity and va the velocity related to the 
amplitude of the superconducting order parameter, the 
chiral symmetry group would be enlarged into the six- 



teen dimensional ?7(4). The SDW insulator corresponds 
to breaking of this symmetry along one particular " direc- 
tion", while the broken generators of J7(4) rotate it to- 
wards one of the other possible broken symmetry states. 
Among these, one can discern three additional types of 
insulators, related to the SDW (and to each other) by chi- 
ral symmetry: the ^'d + ip" and "d + is^' insulating states, 
and the stripe-like one-dimensional charge density waves 
(CDW), in the classification of refs. ^ and ||]. 

In the case of the maximal U{A) chiral symmetry, all 
the above insulating states, and their various combina- 
tions Q|, are equally likely outcomes of the spontaneous 
symmetry breaking. This degeneracy is in reality re- 
moved by the symmetry breaking perturbations to the 
QED3, most prominent of which have already been men- 
tioned above. For example, it was shown in that the 
short range repulsion between quasiparticles favors the 
SDW over the d + ip insulator, and, moreover, enhances 
the SDW order deeper in the insulating state. In this 
paper we study the effects of weak velocity anisotropy, 
higher order derivatives, and short-range repulsive in- 
teraction on the pattern of chiral symmetry breaking in 
the QED3 theory of the phase fiuctuating d-wave su- 
perconductors. All of these perturbations are irrelevant 
at low energies, and we show that their effects on the 
energies of various states with weakly broken chiral sym- 
metry can be ordered according to their engineering di- 
mensionality: velocity anisotropy, being only marginally 
irrelevant Q], |^], provides then the dominant perturba- 
tion. Weak repulsion and the second order derivatives 
are both equally (linearly) irrelevant, but to the first or- 
der it is only the repulsion that affects the energies of the 
insulating states. Fine structure of the chiral manifold of 
insulators is schematically depicted in Fig. ^. 

Our main conclusion is that, to the leading order, it is 
the SDW state that is selected by the weak perturbations 
to the QED3 . Although in real systems none of the above 
perturbations is truly weak, we believe our result pro- 
vides a useful qualitative guide. In particular, it agrees 
with the standard picture of underdoped cuprates, upon 
identification of the SDW insulator as being continuously 
connected with the Mott insulating antiferromagnet near 
half filling. 

In what follows we first briefly review the salient points 
of the QED^ theory of phase fluctuating d-wave super- 
conductor, and set up the convenient notation in terms 
of the eight component Dirac flelds. We then proceed to 
calculate the lowest order splittings of the energies due 
to velocity anisotropy, short-range repulsion, and higher 
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order derivative terms. We end with the summary and 
some concluding remarks on our results. 



II. QEDs 

We start with the finite temperature (T ^ 0) quantum 
mechanical action for the interacting d-wave quasiparti- 
cles, 



S — Sbcs + Si 



(1) 



where 

•^Bcs = 



(2^ 



+ A(fc)cj (fc, W„)cj {-k, -UJn) + h.c. 



with ujnS denoting fermionic Matsubara frequencies, and 
with A(k) as the standard d-wave gap. Su represents the 
short-range repulsion between electrons, 



u 



with U > 0. Shifting the momenta as k = Ki + q,i — 
/, //, and rotating the coordinate frame as in Fig. ^ we 
define the eight-component Dirac field = (^/ , ^/7^) 
where for each node i = I, II, 



(4) 



At r = 0, with qo = Lv, the action may be then written 
as 



d^q 
(2^ 



^(g)ro{z(7o + iM(g)}*((z) 



+ U J ^^r^^2^(^'(r,r)Bc/^'(f,r))^ 

with M{q) = diag{Mj{q^, qy), Mij{qy,qx)), where iMi 
diag(7ii, Hj) is a 4 x 4 matrix defined as 



(5) 



Hi 



i{K, + q) A{K, + q) 



(6) 



Here ^ ~ ^'^Tq, and Tq = diag(7o,7o), with 70 = 

CTl (8) l,7i = -tT2 <8) 0-3,72 = -(T2 ® cri,73 = -0-2 ® Cr2, 

satisfying the CHfford algebra {7^,7,^} = 25 Finally, 
Bu = TqAu, with Au = diag(cr3, 0-3, 0-3, 0-3). 

One may next expand iM{q) around the nodes of the 
superconducting order parameter in powers of q. It is 
convenient to define the following matrices: Mi = —ia^® 




^2) FIG. 1. Four nodes of the d- wave order parameter. Dashed 
line is the putative Fermi surface. 



CT3, M2 = — jT3 ® ai. Using the symmetry property of the 
d-wave gap and the quasiparticle dispersion A(— fc) = 
A(fc) and £,{—k) = ^{k), one finds: 



(3) -S-Bcs = So + Sa + S, 



NL, 



(7) 



So = J drd'r^ir, T){Todr + Fi^, + T2dy}-^{r, r), (8) 
Sa^ f drd^r^if, r){Bi9, + B29J*(r, r), (9) 



where F, = TqA^B, = ToSA^Ai = diag(Mi,M2), 
A2 = diag(M2,Mi),Mi = diag(AFMi, AAM2), M2 = 
diag(AA-M2, XpMi) and Xp ^ vp — 1, Aa = wa — 1 with 
'"F = dq^^{Ki + q)\q=o,VA = dqyA{Ki + g)|q=o- VF and 
va are the two characteristic velocities of the linearized 
spectrum, in units of some fixed velocity c, which we set 
to c = 1. The isotropic theory is recovered in the limit 
vp — UA. We will discuss the form of the Snl that rep- 
resents the higher order derivative terms shortly. Note 
that the matrices ^^'s and SAi's are independent of Fq. 

In the rest of the paper we will consider the linear, 
isotropic, non-interacting Sq to be the maximally sym- 
metric action. It remains invariant under a global chiral 
rotation 



(10) 



with the generators Ji,i — 1 . . A6, forming the ?7(4) al- 
gebra. The explicit form of the generators is given in [0 . 
The connection to the QED3 is finalized by coupling 
the gapless quasiparticles to the fluctuating vortex ex- 
citations, which upon a singular gauge transformation 
may be represented by an effective gauge field in Sq, by 
df_i ^ + ia^. The charge of the gauge field is pro- 
portional to the vortex condensate Such a gauge 
transformation also turns the d-wave quasiparticles into 
neutral spin-1/2 excitations. In the phase coherent, i.e. 



2 



superconducting state, charge vanishes and the gauge 
field is effectively decoupled, leaving the sharp spinon (or 
quasiparticle) excitations behind. When the vortices con- 
dense, on the other hand, superconductivity is lost, and 
simultaneously the spinon "mass" term ~ m^'i' becomes 
dynamically generated Here, m — E(q 0) ^ 0, 
with representing the spinon self-energy due to the 
integraction with the gauge field. Rewriting the mass 
term in terms of the electron operators one recognizes 
that it simply represents a two-dimensional SDW with 
the periodicity determined by the wave vectors 2Ki 

The SDW order, therefore, corresponds to the partic- 
ular form of the Fq matrix we have chosen. By applying 
a unitary U (4) transformation on our Tq one recovers all 
other insulating states that are related to the SDW by 
chiral symmetry. In particular, transforming Fq as 



Fo 



Toe" 



(11) 



with J = 1 (g) 73 produces the d + ip insulator: the 
phase incoherent, translationally symmetric state with 
an additional p-wave pairing between spinous. Similarly, 
J = ^cri(8i[(7i-72 + i7o(7i+72))] yields the d -I- is insu- 
lator, and J — -^0-2 ® (71 — 72) gives the CDW with the 

periodicity set by the wavevector Kj — Kjj. Of course, 
there is a continuum of other chiral rotations that are 
possible, and which lead to various linear combinations 
of the four fundamental states defined above. 

In the rest of the paper we study how the degeneracy 
between the SDW, CDW, d + ip, and d + is states as 
defined above is removed by the leading perturbations to 
Sq: anisotropy, repulsion, and higher order derivatives. 

III. ANISOTROPY 

We set U — in S to study first the effect of weak 
anisotropy (Af, Aa 1) on the energy of the degenerate 
ground states of the chiral symmetry broken, isotropic 
QED^. Since the action S'bcs is symmetric under the 
exchange vp ^ ^a, the energies of various states should 
be invariant under the same exchange. Moreover, if 
vp = VA, both velocities can be rescaled out of the prob- 
lem by an appropriate choice of c. The energy split- 
tings between the states must therefore be proportional 
to (Af — Aa)^, to the lowest order. Indeed, defining the 
three dimensional volume as ^ = {27t)^^ J drd^r, the 
first order correction of the energies per volume due to 
the velocity anisotropy is 

A£;W = l(5i)o = ^tr(S,F,)/o 



V 

Sm^iXp + Aa)/o, 



(12) 
(13) 



independently of the choice of Fq. 
mensionless integral 



Iq is a positive, di- 



/n = 



1 /-^Z" d^3 



(14) 



where x = q/m and cr(x) — I](x)/m is the rescaled self 
energy. A denotes the upper cut-off. ^e'^^ therefore 
provides only an overall energy shift, same for all states. 
To the second order, however, we find 



(2) 



AE 



where 



{{s!)o - {Si)l) 



1 

'w 



1 



-m 



(15) 



^/"^ d^x a\x)x^Xj 
(27r)3 {a^x)+x^)^'^ 

(2^)3 (a2(4+x2)2' 



(16) 
(17) 



Since the rescaled self-energy a{x) in the QED^ falls off 
rapidly for a; 3> 1 Q , we will set A = m in /y , and then 
approximate (t(x) — 1 under the integral. This "constant 
mass" approximation seems not to be appropriate for the 
integral lij^afS, which is without cr(x) in the numerator of 
the integrand. Luckily, however, tT{BiTaBjTp) is inde- 
pendent of the choice of Fq and therefore, the term with 
lij.ap does not contribute to the energy differences, but 
only to the overall shift of the energies. This will be the 
generic situation in all further calculations. We will then 
take A = TO and a{x) = 1 in all the integrals hereafter. 

The lowest (second) order effect of the velocity 
anisotropy is to increase the energy of the CDW relative 
to the SDW, d-|-ip and d+is, which remain degenerate: 



AE 



A, CDW 



AE 
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A,othcr = Am^ivF - VAfi > 0, (18) 

d^x x^ _ 10-37r 
(27r)3 (H-a;2)2 ~ 48tt^ ' 



The energy splitting vanishes when vp = va, as ex- 
pected. 

The fact that d + ip and d + is insulators have the 
same energy in presence of velocity anisotropy can be 
shown to be generally true to any order of the perturba- 
tion theory. To see this, note that d + ip and d + is states 
are represented by To^d+tp = diag(cr2, -0-2, o'2, -0-2), and 
Fo,d+js = diag((T2, (72, (72, f2), respectively. When writ- 
ten in the 2x2 from, the only difference between the 
two is in the signs of some terms, which always may be 
changed by a unitary transformation. Put differently, the 
choice of Fq enters the energy calculation only through 
the combination Bi — ToSAi. Matrices Bi, on the other 
hand, have to appear in even numbers in our calculation, 
as in Eq. (|l5|), otherwise the accompanying integrals will 
be odd in some momentum component and vanish. Since 
Ai 's are block-diagonal, the sign of the block-diagonal el- 
ements in Bi's, then, can not matter: d + ip and d + is 
states remain degenerate to all orders in weak anisotropy. 
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IV. REPULSION 

We now set vp — va — ^ to work out the first finite 
energy contribution of the short-range repulsion Su to 
the degenerate ground states of the isotropic action, ft 
is found that 



AE, 



J = 



(1) _ 



i-{Su)o - m^UJ^ mBu))' - triBfj)], (f9) 



V 

(fx 



(27r)3 f 



4 - TT 

87r2 ■ 



The first term in Eq. vanishes identically for all 

states; the second term also vanishes for the CDW, but 
it increases the energy of the d + ip and the d + is, while 
decreasing the energy of the SDW 



AE, 



(1) 

UA+ip 



ae: 



(1) 

U,d+is 



ae: 



(1) 

C/,SDW 



AE 



(1) 

(7, CDW 



0. 



(20) 
(21) 
(22) 



Note that d + ip and d + is remain degenerate in presence 
of the repulsive interaction as well. 



V. HIGHER-DERIVATIVE TERMS 

We may expand iM{q) beyond the linear terms that 
led to the anisotropic action, Eqs. (||, ||). For instance, 
defining = 1(E) a^, Ma = 1 ® cti , the second-derivative 
term is given by 

Snl = - J dTd^r^jHM^eoid^) + MAA'i{d^)}^i 

+ {I ^ II,x^y} + ..., (23) 

where Cq{xi,X2) = XtXjdxidx^^ixi,X2)\s=o and sim- 
ilarly for Aq. Ellipsis indicates higher than second 
derivative terms. It seems, then, that we need to specify 
the dispersion relation and the gap function to determine 
the effect of the higher-order derivative terms on the en- 
ergies of the degenerate ground states. Interestingly, it 
turns out that up to the the first non-vanishing (second) 
order of perturbation, the qualitative effect of these terms 
is to raise the energy of the SDW, lower that oid+ip and 
d + is, and keep the energy of the CDW unchanged, in- 
dependently of the functional form of the dispersion and 
the gap function. 

For definiteness, we present here the results for the 
tight-binding model on a square lattice, for which the 
dispersion and the d-wave gap are given by 

^(fca, fcb) — — 2t(cos kaCL + cos kbu) — II 
A{ka,kb) = Ao{cos katt — cos ki,a), (24) 

where t is the nearest- neighbor hopping matrix element, 
a is the lattice spacing, /z is the chemical potential, and 



Aq is the amplitude of the d-wave superconducting order 
parameter. The result is then 



A£:l?| , , , = A£;'^^ 



-8m^{ma)^L, (25) 
^JVL,SDW - +&rn^imafL, (26) 

(27) 

Again L is a positive, dimcnsionless integral, given by 



NL.d+ip ^-^NL,d+is 

AE^'' 



A p(2) — n 



L = 



[coikpf r d^ 



1 



4(ma)4 J (27r)3 (l + a;2)2 
X {\vF{coa{max'^) + cos{max^) — 2)]"^ 
+ [vA{cos{max^) — cos(maa;^))]^ } , (28) 

where kp — -^kpa and x^ — -^{Qx iqyj/rn. Assuming 
that near the transition, ma <C 1, to the zeroth order in 
ma yields. 



1 . r d^£ 



1 



L = -{coikpfj (2,)3(i + ,2)2 

X [vl{xl + xlf + v\{2xiX2f] +0{{maf) 
L 
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{coikFY{2vi, + vi;) + 0{{maY), (29) 



where L = (IStt - 46)/1807r2 



VI. SUMMARY. 

Collecting our results together, the shift in energy of a 
given state may be written as 

AE = 8m^ I {vf - va?! + Ou {mU}J^ 



+9nl {ma^L 



(30) 



where each 9a,u,nl = or ±1, depending on the change 
in energy of the state in consideration. In our constant 
mass approximation all three integrals /, J and L are 
mass independent positive constants, so that the Eq. ( |30| ) 
represents an expansion of the energies of the insulating 
states in powers of the dynamically generated mass m. 

For weak chiral symmetry breaking, there is a hierar- 
chy of the perturbation terms, in which each symmetry 
breaking perturbation assumes a place according to the 
degree of its irrelevancy. Velocity anisotropy then, be- 
ing the marginal perturbation at the bare level, provides 
the dominant contribution, while the repulsive interac- 
tion, is the subdominant one. Second order derivatives, 
although equally irrelevant as the repulsion, contribute 
to the energies of the insulating states only to the second 
order, and therefore are the least important perturba- 
tion. Breaking of the degeneracy between the insulators 
due to each perturbation is schematically given on Fig. ||. 
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SDW,dip,dis,CDW 



Bare Theory 




Anisotropy 



Higher-Order 
Derivatives 



FIG. 2. Energies of the four types of insulators in the un- 
perturbed (bare) action, with the velocity anisotropy, repul- 
sive interaction, and the non-linear (higher order derivative) 
terms. 



Our conclusion is that the degeneracy of the chiral 
manifold is broken in favor of the SDW, which is the 
lowest energy state when the weak perturbations to the 
maximally chirally symmetric action are taken into ac- 
count. The translationally symmetric d + ip and d + is 
insulators remain degenerate, as one would expect, since 
the translational symmetry remains intact even in pres- 
ence of all three perturbations. Assuming that d-wave 
superconductor-insulator transition is continuous, or pos- 
sibly weakly first order, implies then that the insulating 
state is the SDW. If the transition is strongly first order, 
chiral mass m increases and it is conceivable that there 
could be some level crossings in our fine structure of the 
chiral manifold. 
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